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We develop a relativistic variational model for a nematic liquid crystal interacting with an electro- 
magnetic field. The constitutive relation for a general anisotropic uniaxial diamagnetic and dielectric 
medium is analyzed. We discuss light wave propagation in this moving uniaxial medium, for which 
the corresponding optical metrics are identified explicitly. A Lagrangian for the coupled system of 
a nematic liquid crystal and the electromagnetic field is constructed, from which a complete set of 
equations of motion for the system is derived. The canonical energy-momentum and spin tensors are 
systematically obtained. We compare our results with those within the non-relativistic models. As 
an application of our general formalism, we discuss the so-called Abraham-Minkowski controversy 
on the momentum of light in a medium. 

PACS numbers: 03.50.Dc; 47.75. +f; 83.80.Xz; 47.10.ab; 03.50.-z 



I. INTRODUCTION 



Liquid crystals provide an interesting example of a 
subject where the fundamental and applied sciences are 
deeply related. After their first experimental discovery 
more than 120 years ago, great number of substances 
(natural and synthesized) with the properties of liquid 
crystals are known, which have many important practi- 
cal applications for modern technology. Good overviews 
and introduction to this subject can be found, for ex- 
ample, in [1-8]. In our study, we deal with nematic liq- 
uid crystals (with a possible generalization to cholesteric 
crystals) which fall into a particular class of media with 
microstructure. In classical continuum mechanics, a ma- 
terial medium consists of structureless points. In the 
early 20th century, the Cosserat brothers [9] proposed a 
generalization of this simple picture, in which the ma- 
terial body or fluid is formed by particles whose micro- 
scopic properties contribute to the macroscopic dynamics 
of the medium. These more complex continuous mechan- 
ical models are known under different names such as the 
theories of multipolar, micromorphic or oriented media 
[10-12]. An important particular case of continua with 
microstructure is represented by the spinning fluids [li^- 
15]. 

As a medium with anisotropic electromagnetic (opti- 
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cal) properties, a nematic liquid crystal is another partic- 
ular case of a medium with microstructure. Just like the 
spinning fluid, which is characterized by elements with 
an "internal" degree of freedom associated with spin, the 
liquid crystal is a medium of "stretched" particles whose 
orientational motion is described by an additional hydro- 
dynamic variable. Mathematically, this additional degree 
of freedom is a unit vector field n — n{x,t), ri? = 1, 
which is called director. The director is a microscopic 
variable that is assigned to every material point of the 
medium. For the cholesteric liquid crystals, in addition, 
the chirality (handedness) property is assigned to the ma- 
terial points. 

We construct here a complete relativistic Lagrangian 
theory of nematic liquid crystals interacting with the 
electromagnetic field. The non-relativistic variational 
models were developed previously in [l(j, IT], see also 
[7] . This variational approach is convenient for the study 
of the full nonlinear dynamics of a liquid crystal, namely 
the equations of motion and the conservation laws. Rel- 
ativistic fluid models are working tools in various fields 
of research such as high-energy plasma astrophysics and 
nuclear physics (where non-ideal fluids are extremely suc- 
cessfully applied to the description of heavy ion reactions) 
[18, 19]. Also in cosmology, hydrodynamical descriptions 
of matter are standard both for the early and for the later 
stages of the evolution of the Universe. Our derivations 
make use of the earlier studies in which the relativis- 
tic Lagrangian theories were developed for ideal fluids 
without [20-22] and with microstructure [2.'^)]. Of special 
interest are the models [24, 25] of relativistic spinning 
fluids. Since we use the Lagrangian formalism it is non- 



trivial to take into account dissipativc effects. In this 
work we neglect dissipation in the motion of the rela- 
tivistic liquid crystal and therefore our model has to be 
understood as a first step towards a more realistic theory. 
On the other hand, the description of light propagation in 
this medium is not affected by the inclusion of viscosity. 
Additionally, even in the case in which the medium is a 
"rigid" non-dissipative anisotropic (bircfringent) crystal, 
a model in terms of a "liquid" (i.e. a fluid) is needed to 
consistently describe its response to the electromagnetic 
field, since the concept of rigid bodies is incompatible 
with a relativistic description. 

The general framework of the latter works has been 
used recently for the investigation of the problem of the 
energy and momentum of the electromagnetic field in 
moving media [2G, 27]. However, this study has been 
restricted to isotropic media only. Here we further de- 
velop a relativistic Lagrangian fiuid model and apply the 
extended theory to the description of a nematic liquid 
crystal and its interaction with the electromagnetic field. 
This is an example of light interacting with a dynamical 
anisotropic medium, which can be used to gain a deeper 
insight in the long standing Abraham-Minkowski contro- 
versy [28-34]. 

As a starting point we use the nonrelativistic model 
[16, 17]. The important difference is however that we use 
the Euler picture and not the Lagrange one as in [16, 17]. 
The Euler approach seems to be more physically attrac- 
tive since it provides the field-theoretical treatment for 
both the electromagnetic field and the material variables 
that describe the dynamics of the liquid crystal. 

Our discussion here is confined to the flat Minkowski 
spacetime. In other words, we do not consider the gravi- 
tational effects encoded in a nontrivial Riemannian met- 
ric. However, the generalization of the relativistic La- 
grangian model of liquid crystals to curved spacetimes is 
fairly straightforward. A possible physically important 
application of such a generalization would be a reanaly- 
sis of the structure and defect formation (walls, strings, 
textures) in cosmology, relying on the analogy between 
cosmological defects and defects in liquid crystals [35, 36]. 

Our notation follows [2(i] and the book [37]. In partic- 
ular, the indices from the middle of the Latin alphabet 
i,j,k,... = 0,1,2,3 label the 4-dimensional spacetime 
components, the Latin indices from the beginning of the 
alphabet a, 6, c, ... = 1,2,3 refer to the 3-dimensional 
spatial objects and operations (the 3-vectors are also 
displayed in boldface). The Minkowski metric is de- 
fined as gij := diag(c^, —1,-1,-1) and the 3-dimensional 
Levi-Civita symbol in the rest frame is chosen so that 
gi23 ._ 2. The determinant is denoted as usual by 
g = det{ gij); thus in Minkowski spacetime y^—g — c. In 
the 4-dimensional framework spatial components of ten- 
sor must be raised or lowered by gab = —Sab, but when 
we are working only with 3-dimensional tensors, we use 
the convention of using just the Euclidean metric Sab to 
raise and lower spatial indices. 



II. CONSTITUTIVE RELATIONS OF A LIQUID 
CRYSTAL AT REST 



In Maxwell's theory [37] the electromagnetic field is 
described by the electromagnetic field strength Fij = 
{E, B) and the electromagnetic excitation H^^ = 
{D,H). Additionally, in order to obtain a predictive 
theory, one needs to specify the constitutive relations 
H^^ — H^^ [Fi^i) for the specific medium under consid- 
eration. 

A liquid crystal is a medium with uniaxial anisotropic 
properties. The constitutive relations for a medium with 
electric and magnetic properties of this kind can be ex- 
pressed, when the medium is at rest, by 



7-ja ^ °a6 



eo^" Eb, 
Ha^t^o'kbB". 



(2.1) 
(2.2) 
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Here £" is the relative permittivity tensor and M^^ is the 

o 

inverse of the relative permeability tensor fJ- . The fact 
that all quantities are considered in the frame where the 
medium is at rest is indicated with the symbol °. 

Alternatively to (2.1) and (2.2), it is also useful to write 
the constitutive relations in terms of the polarization and 
magnetization fields: 



P'^ = eo(|"'' - 5'''')Eb, 

Ma^ll^\5ab~Kb)B'- 



(2.3) 
(2.4) 



In the case when the medium has the same optical 
axis n for the electric and magnetic anisotropy, it is con- 

o . ° — 1 

venient to decompose e and M^,^ in terms of its eigen- 
vectors. By choosing these vectors as n, ni, and n2, with 
two eigenvalues equal and one different, we have 



£°'' = eiin'^n'' + i 



/ a b , a b\ 



M, 



ab 



/in ^naUb + Hj^(nlnl + ri^ng) 



(2.5) 
(2.6) 



Here e± and eu are the relative permittivities, perpen- 
dicular and parallel to the optical axis vector n (with 
Cartesian components n""); ii±_ and /xy are the corre- 
sponding perpendicular and parallel relative permeability 

o , ° — 1 

functions. Since the tensors e" and M„{, are symmetric 
and real for a non-dissipative medium, the eigenvectors 
satisfy the closure relation 



a b , a b , a b 

n n + niUi + 7^2712 



gab 



(2.7) 



Inserting (2.7) in (2.5) and (2.6), we obtain a simpler 
expression for the dielectric and magnetic tensors when 
the medium is at rest, 



iab 



cafc 



Ae n°n' 



a^b 



M_L^ Sab + A/i ^Ua 



rib- 



(2.8) 
(2.9) 



The dielectric and magnetic anisotropics are defined, re- 
spectively, by 



Ae: 
An~^ := ^li 



1 



Mj 



(2.10) 
(2.11) 



The latter quantity should not be misunderstood as the 
inverse of the difference /X||—/x_L, that is, A/i~^ ^ (A/x)~^. 
Strictly speaking, one should write (2.11) as A{^~^), but 
we omit the parentheses to simplify the formulas. 



III. NON-RELATIVISTIC LAGRANGIAN FOR 
A NEMATIC LIQUID CRYSTAL 

The non-relativistic liquid crystal theory is a well es- 
tablished subject; see, for instance, [1, 3, 8, 38-44]. Nev- 
ertheless, the Lagrangian approach for the study of the 
dynamics of this medium with microstructure was devel- 
oped only recently in [16, 17] (although the variational 
methods were used in [7] for the analysis of the equilib- 
rium problems for the liquid crystals). 

According to [16, 17], the non-relativistic kinetic en- 
ergy density of a liquid crystal reads 



K. := p, 



V 



PmJ- 



(3.1) 



Here pn^{x,t) is the mass density of the liquid crystal, 
v{x,t) its velocity field, J^ the geometric moment of iner- 
tia of a fluid element (dimensionless) , and lj the angular 
velocity of the director, which is defined as 

uj -.^ n X n, (3.2) 

with n{x,t) the director field of the liquid crystal and 

dn 



dt 



{v ■ V)n, 



(3.3) 



the convective derivative of the director. 

The potential energy density is represented by the free 
energy J-, which is usually taken as the thermodynamic 
potential in the theory of liquid crystals. Following [1, 8, 
16, 17], we express the free energy as 



J^ — J^o + J^d + J'c + J^n 



(3.4) 



where J"o = ^o{Pin,T) is the internal free energy which 
describes the hydrodynamic portion of J" and depends on 
the density p^ and the temperature T. The pressure in 
the medium is introduced as p := p^^^{dJ-'o/dpin)T- The 
internal dynamics of the director field is described by 
the Frank deformation potential J-d which is defined in 
[1], for the simpler liquid crystal with group symmetry 
{oo/ mm), by 



JFd = :^Ki{V -n 



)2 + ]-K2{n • V X n)2 



-K:i{n X V X nf . 



(3.5) 



The three parameters Ki, K2 and K^ are known as 
Frank's elastic constants (elastic moduli), which are all 
independent from each other and also positive. One usu- 
ally calls Ki splay, K2 twist, and K^ bend constants. 
The so-called saddle-splay boundary term is omitted, 
since it is a total derivative that does not contribute 
to the equations of motion. For a typical nematic crys- 
tal, one has Ki = 2.3 x IQ-^^N, K2 = 1.5 x IQ-^^N, 
Ki ==4.8x 10-12 N (see [1, 7]). 

The cholcstcric liquid crystals arc characterized by an 
additional modulus A'o and a constitutive constant r, 
related to the chirality of the medium. As a result, the 
Frank potential (3.5) is modified to 

Fd = A'or(n -V y.n + T) + -A'i(V • nf 

+ -K2{n -V xn + rf + -K^{n x V x nf. 

(3.6) 

We restrict ourselves to the case of the nematic crystals 
with T — 0, although the generalization to the cholesteric 
crystals is straightforward. 

The interaction free energy of the liquid crystal with 
an electric field E is represented by F^- Generally, con- 
trolled in the field E, the electric free energy of the sys- 
tem [45] reads J^o = — J P ■ dE. Using (2.3) and the 
expression for the permittivity tensor in the comoving 
frame (2.8), we explicitly obtain 

Fe = -^eo(ei - 1)^' - IsoAein ■ Ef. (3.7) 

Analogously [3, 46], we have J"m = — J M ■ dB for the 
magnetic free energy controlled in the field B, which 
yields 



Tr, 



J-{^l-pl^)B^ + -^Ap-\n-Bf. (3., 
2p,Q 2pa 



Then, the nonrelativistic Lagrangian of the nematic 
liquid crystal is constructed as the difference £„!■ = IC—T 
of the kinetic energy density /C and the "potential" free 
energy density T. Accordingly, we have 



v^ U}^ 1 

Air = Pm-^ + PmJ'-^ - J^o{Pm,T) - -Ai(9a"-") 



a\2 



^K2ie''''^nadbn,f - ^K^ieabcn^e^^'Ban^f 



+ -eof±-E' + -£oAe(n • Ef 

--}-pl^B'-^Ap-\n-Bf, 
2^0 2/io 



(3.9) 



where we added the energy density of the pure electro- 
magnetic field eo-E'^/2 — B'^ /2po. This is necessary to 
describe the electromagnetic field as a dynamical part 
of the system and to guarantee the correct limit of the 
energy density in free space. 



The total non-relativistic liquid crystal Lagrangian 
(3.9) can be conveniently recast into the sum of the mat- 
ter part £™ (pm, T, u°, n°, dbn"-) and the electromagnetic 
part£°-(n^£;„B'^), 



n nm. i pern 



(3.10) 



where 



- ^i^2(e'""=n,47i,)2 - lA'3(ea6c«'e^^^adne)^ 



(3.11) 



i^T = leoe^E^ + ^eoAein ■ Ef 



- ^Mo^Vl'B' - \^^^'^^i-\n ■ Bf. (3.12) 

The non-relativistic variational theory based on the La- 
grangian (3.9) was developed in [IG, 17]. It used the La- 
grange approach and was formulated in terms of {E, H) 
instead of the fields {E^ B). Dynamics of the total system 
is determined by the action integral /m- = J dtcPxCm-. 



IV. RELATIVISTIC LIQUID CRYSTAL 
LAGRANGIAN 



The formal theory of liquid crystals has been studied 
only in the non-relativistic domain, as can be seen in [1~ 
8, 16, 17, 38-41]. Here we develop a truly relativistic 
model for these systems, generalizing the non-relativistic 
3-dimensional objects and operations to the correspond- 
ing 4-dimensional notions. In contrast to [16, 17], we 
work in the Euler approach which is more convenient for 
field-theoretical applications. 

First, we notice that the 3-velocity field of the liquid 
crystal v is proportional to the spatial part of the 4- 
velocity m* = (7,71;), which by definition always satisfies 
the condition 



>0. 



(4.1) 



Thus, u 
1 



1 



a timelike 4- vector field. Here 7 = 
v^/c^ is the usual Lorentz factor. Analogously, 
we can define the director ^-vector TV' as the relativis- 
tic covariant generalization of the director n" . When the 
medium is at rest. A''* should reduce to n"; that is. 



7V' = (0,n). 



(4.2) 



We now recall that n, by definition, has a unit length. 
This together with (4.2) imposes the scalar condition 



N'N, = -1 < 0, 



(4.3) 



which should be fulfilled in all reference frames. In other 

o ■ 

words, N"^ is a spacelike 4-vector. In the rest frame, u^ = 



(1,0), and together with (4.2), we have UiN"^ = 0. Since 
this is a scalar condition, it must be valid in all reference 
frames as well: 



N'ui = 0. 



(4.4) 



In addition, it is necessary to define a relativistic gen- 
eralization of the 3-dimensional Levi-Civita symbol in or- 
der to consistently express the "cross products" in (3.9). 
Let us introduce 



^ijk • Vijkl ' 



(4.5) 



with the 4-dimensional Levi-Civita tensor defined such 
that 770123 := V~9 — c and thus in the rest frame its spa- 
tial components reduce to the usual 3-dimensional Levi- 
Civita symbol e°^'^, with e^'^^ = 1. Using this object, we 
immediately define the angular 4- velocity of the director 
by 



w 



e'^'^N^Nk 



(4.6) 



where the convective "time" derivative is naturally 



N' 



u^djN' 



(4.7) 



The relativistic variational theory of an ideal fiuid with 
structureless material elements is a well developed sub- 
ject [20-23, 26]. The generalization to the ideal fluid 
with classical spin (modeled after the Dirac particles) 
was done in [24, 25]; see also the references therein. 

A liquid crystal medium represents another example of 
a fluid with microstructure, represented by the director 
field attached to each element of the fluid. Here we de- 
velop a relativistic variational model for this system by 
combining the variational model described in [26] (for the 
kinetic translation energy and the internal energy), with 
the four-dimensional generalization of the liquid crystal 
elastic terms in (3.11): 

+ ^Ks {e,jkN^e^'--diNnf + Ao(u'm, - c^) 



+ ki{N^N, + l)+kzu^N,. 



(4.8) 



Here v is the particle number density of the liquid crys- 
tal, J is the moment of inertia of one element [related to 
J^ in (3.1) by means of J Pm = Jv], p{v, s) is the internal 
energy density of the relativistic fluid, s is the entropy 
density, X is the identity (Lin) coordinate, and A/, with 
/ = 0, . . . , 5 are Lagrange multipliers. By imposing these 
A's we ensure the fulfillment of the conditions (4.1), (4.3), 
and (4.4) throughout all the dynamics of the nematic liq- 
uid crystal, in addition to the particle number continuity 
equation. 



d^iuu')=0, 



(4.9) 



and the conservation of entropy and identity of particles 
along each streamline of the fluid: 



u'^diS = 0, 
u'd^X = 0. 



(4.10) 
(4.11) 



A different sign of the K3 term in (4.8), as compared 
to (3.11), is explained by the fact that the 4-vector 
eijkN^e^''^diNn is spacelike, hence the square of its 4- 
length is negative. Notice that, as is usual in relativistic 
fluid models, we choose s and v as the indepedent ther- 
modynamic quantities, instead of T and pm as in the 
non-rclativistic case. This choice is purely convencional, 
since all the other thermodynamic quantities can be de- 
rived from them. 

The dynamics of the relativistic system is governed 
by the action / = (1/c) J y/—gd'^xC'^. One may won- 
der how the nonrelativistic translational Lagrangian can 
be recovered from the relativistic Lagrangian. As a first 
step, we represent the internal energy density as the sum 
P = PmC^ + ^0 of the "rest-mass" density and the hy- 
drodynamic energy density. Consider now an arbitrary 
volume element which reads \/—g d^x — c dVb dr in the 
comoving reference frame with the 3-volume dV^ and the 
proper time r. The next step is to notice that the rest 
mass of a fluid's element drriQ = Pm^Vo is the same in all 
frames and the invariant volume element dVo dr = dVdt 
[with dr = ^1 - v'^/c^dt w (1 - v^/2c^)dt] in the ref- 
erence frame where the fluid's element dV has veloc- 
ity V. As a result, in the nonrelativistic limit we in- 
deed recover the translational part of the Lagrangian: 
(1/c) / y/^d'^x {-p) « / dtdV IpW/2 -To). A more 
detailed discussion can be found in [47], for example. 



V. LAGRANGIAN FOR THE 

ELECTROMAGNETIC FIELD INTERACTING 

WITH THE LIQUID CRYSTAL MEDIUM 

In order to describe the interaction of the electromag- 
netic field with matter in an explicitly covariant manner, 
we make use of the standard electromagnetic Lagrangian 
given by 



£"" = - -W^F, 



(5.1) 



It yields the macroscopic Maxwell equations as Euler- 
Lagrange equations without sources (we assume that the 
fluid elements are not electrically charged): 



djH'^ = 0. 



(5.2) 



Here the electromagnetic strength tensor Fij is expressed 
in terms of the electromagnetic 4-potential Ai, as usual 
by F,, 
magnetic excitation tensor, [26, 37]. 

The constitutive relations for any linear, non- 
dissipative and non-dispersive medium can be expressed 



diAj — djAi and W^ is the covariant electro- 



in general covariant form by 

1 



H-i = -x'^'^'Fu, 



(5.3) 



where x*^*"'' is the so-called constitutive tensor with the 
symmetries 

ijkl _ _ jikl _ _ ijlk _ klij /g ^\ 

Therefore, it has 21 independent components, in general. 
If we insert (5.3) into (5.1), we obtain 



C''''^-lx'''''F,,Fki. 



(5.5) 



Historically, the general constitutive relation (5.3) was 
first formulated by Bateman [48], Tamm [49-51], and 
later in the modern notation by Post [52]. 

As we commented in Sec. II, the nematic liquid crys- 
tal is an example of a uniaxial dielectric and diamag- 
netic anisotropic medium and therefore we need to find 
an explicit covariant expression for the constitutive ten- 
sor x*-**^' for a medium of this kind. 



A. Constitutive tensor in the comoving frame 

The components of the covariant constitutive relation 
(5.3), must reproduce the expressions (2.8) and (2.9), in 

O 7 

the rest frame of the medium. Therefore, given e" and 

° _i . 

P^i^ in terms of the director n and the eigenvalues e±, 

— 1—1 ° ■ kJ 

e||, /.j^ , and fi,, , the non-vanishing components of x*^ 
must explicitly read 



l^OabO 



'ibcd 



Ae 



n n , 



(5.6) 



abf cdg,,-l 



fg 



- Afi-^ {5" 



■-1 ''^''''n''n'' 



S" 



d„b. 



+5'"^n''n''' - (5'"=7i°7i'^) . (5.7) 



It is worthwhile to notice that this constitutive tensor 
characterizes all non-magnetoelectric anisotropic media 
with dielectric and diamagnetic uniaxial properties in the 
direction of the optical axis n. In the special case of 
nematic liquid crystals, the optical axis coincides with 
the director field. 



B. Dispersion relations and factorization of the 
Fresnel equation 

If wc look for wave solutions to the macroscopic 
Maxwell's equations (5.2) in a source-free and homoge- 
neous medium (described by x'"''"' ) , then the general dis- 
persion relation is determined in covariant form by the 
fourth order Fresnel equation for the 4- wave covector fc^: 



g'''''hkjkkki = o. 



(5.8) 



Here tj*-''^' is the Tamm-Rubilar tensor [53], given by 
1 



t'ijkl 



mnr{i j\ps\k l)qtu 



4!c2 



(5.9) 



We can use the non- vanishing components of the con- 
stitutive tensor (5.6) and (5.7) in (5.9) and (5.8) and 
thereby verify the factorization of the 4-rth order Fresnel 
wave surface into a product of two hght cones, deter- 
mined by two optical metrics in the rest frame of the 
medium: 

t^'^'hkjkkh = ill' hk,){^^ kkh) = 0. (5.10) 

Here the hght cones read exphcitly 

%' hkj = n2fc^-aofc^ + (ao -!)("• fc)^ (5.11) 
Tihkj = n2fc2-a,„fc2 + (a,„-l)(n-fc)2,(5.12) 

with the refractive index 



A*J 



(5.13) 



The parameters olc and am quantify the proportion of 
dielectric and diamagnetic uniaxial anisotropy in the 
medium, and read 



ei 




^1 




am := 




e|| 




A^ll 



(5.14) 



One of the reduced quadratic Fresnel dispersion rela- 

o . . 

tions, 1]^ kikj = 0, implies that if fc is parallel to n, then 
n^fcg — fc^ = 0, which means that in this case light prop- 
agates with the expected effective refraction index of the 
ordinary ray: n. On the other hand, if k is orthogonal to 
n then the dispersion relation reduces to n^kQ — ack"^ = 0. 
This means that light propagates with an effective refrac- 
tion index rio, given by n^ = n^/ae = ey/ij.. We may call 
this the dielectric extraordinary ray. Similarly, the sec- 

o . . 

ond Fresnel dispersion relation, "f\l^kikj = 0, leads to a 
normal ordinary ray refraction index n for waves with 
wave vector fe parallel to the optical axis n. For fe _L n, 
it implies that light propagates with a refraction index 
?im , with ?i^j = n^/a-ca = £J-M||i corresponding to the 
diamagnetic extraordinary ray. 



C. Covariant description 

Since both (5.11) and (5.12) must be covariant equa- 
tions, it is not difficult to derive the general expressions 
of the two optical metrics in a frame where the medium 
moves with an arbitrary 4- velocity u' ~ (7,71;): 



-fl^ = ao g' + 



{n? - ac) 



u'u^ + {a^-l)N'N^, (5.15) 



7^ =am5'^ + 



(n^ - am) 



u'u^ + {a^n - l)N'N^ .(5.16) 



One can then verify that the covariant form of the consti- 
tutive tensor, which takes into account all the necessary 
symmetries and which reduces to (5.6) and (5.7) in the 
rest frame of the medium, is given by 



X 



ijkl 



1 

Mo 



M 



-' + A^-i) (g'^g-'"' - g^'g^'^) 



+ - 



MoC 



1 



(n^ - 1) 



/^-L 



A^-1 



(g'''u^u^ - g'^u^u'' 



+g^^u'u'' - g^'^u'u^) 



Mo 

+g3^N'N^ - g'^N'N^) 



i:.(Ae + A^-i) 



I'u^'N^N^ ~ u'u^N^N'' 



Moc" 



+u^u^N'N'' - u^u^N'N^) . (5.17) 



Generalizing the result obtained by Balakin and Zim- 
dahl in [54] , the constitutive tensor can also be expressed 
in terms of the two optical metrics (5.15) and (5.16), by 



ijkl 



1 



Otc l^OlJ-± 



i^/ik^,jl ^M^A^\ 
(7o 7o - 7c 7o J 



1 



(am - ac)/io/i-L 



(A7*''A7-''' - A7''A7-'''') 



(5.18) 
where A7'-' is the difference of the optical metrics, 

A7'^:=7:^-7S- (5-19) 



D. The projector vr'^ and the isotropic limit 

Let us define the projector tt'-' to the 2-dimensional 
space orthogonal to A^* and u*. 



tt'^ ■.= g'^ -\u'u^ +N'N^. 



(5.20) 



It obviously has the properties n^jTr^k = Tr'fe and 
det(7r'-') = 0. With this object we can write the opti- 
cal metrics in a more compact form: 



7^^' = 7'^' + («c - 1)^*^ 
7^ = 7*-'' + ("m - l)7r^^ 



(5.21) 
(5.22) 



where 7'-' is the usual optical Gordon metric of an 
isotropic medium [26, 55], 



Y^ ^ g'^ + 



in' - 1) 



(5.23) 



Also the difference of optical metrics (5.19) can be ex- 
pressed in terms of tt*-' , 



A7'^' = (ae - am)7r*- 



(5.24) 



Therefore, the constitutive tensor in (5.18) can be recast 
into the form 



ijkl 



1 



^^o^^± 



1 



ik jl 



(7c 7, 



7:'7^') 



-(ttc - am) (TT^'^Tr^'' - tt^'tt^'^)] . (5.25) 

From (5.25), we can easily check that our x^-''^^ reduces to 
the well-known expression for the isotropic case in terms 
of the Gordon metric (5.23), when 

e|| ^ £_L = e ^ ac^l, (5.26) 

^J.^ ^ H± = fi 4^ ttni -^ 1. (5.27) 

Then both optical metrics in (5.21) and (5.22) become 
the single Gordon metric, and the constitutive tensor re- 
duces to 

X''"' ^ Xifo' := — (7'S^'' - 7''7^'') . (5.28) 



Then the electromagnetic Lagrangian (5.1) can be alter- 
natively written in a "3-D like" form. 



>C"" = I {eoe'^8^8, A^o V.^'S'S^") 



(5.33) 



where e^^ is the 4-permittivity tensor and /i^ ■ the inverse 
of the 4-permeability tensor, given by 

e'^ := - e^cf^ + /^sN'N^ (5.34) 

Mi,-' := - Mi'ff.j + ^ti~^N,N,. (5.35) 

Notice that the spatial components of these expressions 
consistently reduce to (2.8) and (2.9). Inserting (5.34) 
and (5.35) into (5.33), we can separate the contribu- 
tions of the electromagnetic Lagrangian in "isotropic" 
and "anisotropic" parts, that is. 



C 



cm /'ISO 



£" 



(5.36) 



E. Explicit expression for the electromagnetic 
Lagrangian 

With the help of the constitutive tensor (5.17), we can 
use (5.3) to obtain an explicit expression for the electro- 
magnetic excitation W^ : 

H^^ = — (^tT^ + A/i-i) F'^' + — A^-i i^['=„7V'l7V" 
Mo Mo 



with 



Moc 
2 



(Ae + Am"^) A^I'^w'l FpgNPu". (5.29) 



From here on, we assume that the permittivity and 
permeability can be functions of the particle density. 
Mi = M±(^)i £± ~ £±i^): a-nd A^^^ = Am^^(i^), in gen- 
eral. Then, from (5.1) and (5.29), it is straightforward 
to compute an explicit expression for the electromagnetic 
Lagrangian £'''" = C"''^{v,u^,N^,Fij), which reads 

£<=- = - -^(mI' + Am-^)F,,F'^" - -^Am-1 (FkiN'f 
4mo 2mo 



1 



l\2 



2fiQC^ 



(£± - Ml' - Am"') {Fkiu') 



+ -^ (Ae + Am-^) {Fp^N^u'^f. 



(5.30) 



For completeness, we also give an alternative deriva- 
tion of (5.30) based directly on the non-relativistic La- 
grangian (3.12). In the four-dimensional relativistic 
framework, the electric E and magnetic B fields are sub- 
stituted with the 4- vectors of electric field £i and mag- 
netic field i3', defined as 



f — p. ,,j 



1 



g. := ^rt^^iF.kUi. 



(5.31) 
(5.32) 



£"° = -\ [eoe^E'£^ - ii^\il^B'B,] 



(5.37) 



/:^'" = - [e„/\e{£,N''f - Mo 'AM-'(S,iV')'] ■ (5.38) 

As we see, the structure of (5.36), (5.37), and (5.38) 
follows exactly the structure of the non-relativistic La- 
grangian (3.12). By inserting (5.31) and (5.32) into 
(5.37), we can write the isotropic electromagnetic La- 
grangian D^° in terms of F^^ : 



D' 



1 



-g'''g^'F,,Fki 



4moM 

2moM±c- 

4MoMi 



^g'^u^u'F.jFki (5.39) 



1 



zj/c/ 






(5.40) 
(5.41) 



where xHo i^ ^^^ constitutive tensor for the isotropic 
case given in (5.28) and written in terms of the Gordon 
metric (5.23). In this case n^ = e±fi± as expected. 

Now we do the same for the anisotropic part of the 
electromagnetic Lagrangian (5.38). First, using the defi- 
nitions (5.31) and (5.32), we explicitly have the following 
expressions: 

{Nm^f ^ \f^^F,, + {F,,N^f ^{F,yr 

- l(F,iVV-'y, (5.42) 

{N'£,f = {F,,N'u^f. (5.43) 

Then, using (5.42) and (5.43) in (5.38), we derive the 
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explicit expression for the anisotropic Lagrangian, 

^ ^^(Ae + A^i-')N'N''u^u'F,,Fkl (5.44) 



2fj,oc- 
with 






(5.45) 



Moc- 



+ A^o ^ ^M"^ {g'^N'N^ ~ g'^ N^ N^ 

+ g^^N'N^ - g^^N'N^) 
1 



Moc' 



(Ae + A/i-i) [u^vf'N^N^ - u'u^Nm'' 
+ u^u^N'N^ ~ u^u^N^N^) . (5.46) 

Finally, the total constitutive tensor -x_^^^^ in (5.17), 
is formed by adding (5.18) and (5.28), as can be easily 
checked, 

X'^'' =x:^'(£±,A^±)+Xai^'(A£,AM-^). (5.47) 

The resulting constitutive tensor is identical to (5.17). 



The Lagrangian part containing the constraints reads 
C^ := Ao(w'w. - c^) - vu'diKi + Aau^a^s 

+ K:iu'd^X + KiiN'^N, + 1) + K^u'N,. (6.5) 

A. Kinetic term 

To begin with, we notice that making use of the defini- 
tion (4.6), we can recast the square of the angular velocity 
into 

2 



lo'lo, = P'j NiN^ ^ N'N, 



^{^^N^y^ (6-6) 



where P^ j is the usual projector operator, perpendicular 
to the 4-velocity field: 



P',:^S]-^u^u,. 



(6.7) 



Therefore, the kinetic Lagrangian (6.2) depends only 
on i^,u'-, and the derivatives of the director field dkN^. 
The derivative with respect to the latter reads 

ggJ^-^J-U^P^^N,. (6.8) 

As a result, we find the variational derivatives of the ki- 
netic Lagrangian with respect to its arguments 



(6.9) 



Oz/ oi' 2 

« = a7 = •'" (-^''ft" + "•"'/=') "j. (6-10) 



VI. VARIATION OF THE MATTER 
LAGRANGIAN 



B. Potential term 



With a well defined Lagrangian (4.8) for our relativistic 
model of a nematic liquid crystal, we can now derive the 
field equations of the system. For this task, we write the 
matter Lagrangian £™ = C'^{v,u\N\djN\s,X,Ai), 
/ = 0, 1, 2, 3, 4, 5, as a sum of three terms. 



C'^ = C^ + CP + C^, 



where 



Jiy (jj^LOi 



is the kinetic Lagrangian and 

£P:=-p(z.,,s)-V 



.1) 



(6.2) 



(6.3) 



is the Lagrangian of the potential energy, with the Frank 
deformation potential given by 



-ii^3(eyfe7V-'e''"<9,iV„)'. 



(6.4) 



The potential Lagrangian (6.3) depends on z/, s, m*, A^*, 
and the derivatives of the director field dkN^ . The deriva- 
tive with respect to the latter reads 

dCP ^ dV 
ddkN' ddkN^' ^ ' ' 

We straightforwardly compute the variational derivatives 
SO' _dCP _ dp _ p + p 
5v dv 



du V 

SCP _dCP _ dp _ 

6s ds ds 

SCP _ a£P _ dV 

5CP_ _ SV_ _ dV_ 



( dV 



(6.13) 
(6.14) 
(6.15) 
(6.16) 



In (6.13) and (6.14) we used the thermodynamic (Gibbs) 
law 



Tds = d 



{^)Ml 



3.17) 



C. Constraint term 

Variation with respect to the Lagrange multipliers A/ 
yields the complete set of constraints: 



U Ui = c , 
^^{VU')=Q, 

u'diS = 0, 
u'diX = 0, 

N'N, = - 1, 
N'u, ^ 0. 



(6.18) 
(6.19) 
(6.20) 
(6.21) 
(6.22) 
(6.23) 



Additionally, the variations of the constraint Lagrangian 
with respect to the field variables read 

SC^ dC 



5v 


dv 


(6.24) 


Ss 


\dd,s) 


(6.25) 


6X 




(6.26) 


6C^ 


= -TT— == 2AoWi - i^diAi + K2diS 
ou'- 






+ A35,X + A5A^„ 


(6.27) 


6N' 




(6.28) 



D. Field equations 

We are now in a position to write the field equa- 
tions, collecting the variations of the three terms in (6.1). 
We can distinguish three groups of equations. The first 
group describes the constraints (6.18)-(6.23). The sec- 
ond group concerns the variables s, X which do not enter 
the Lagrangian of the electromagnetic field. This yields 






(6.29) 
(6.30) 



The third group of equations is obtained from the vari- 
ations with respect to the essential field variables ly, m*, 
and A^\ We have explicitly 



3.31) 



—— =- -uj' ^-u'a.Ai, 

ou 2 u 

+ 2AoUi - udiAi + A2dis + A^diX + AsN^, 

(6.32) 

(6.33) 



6N' 



=4 (^Jv u'^P^.Nj)- ^+ 2A^N, + Agw, 



Contracting (6.33) with u^ and A''*, we find the Lagrange 
multipliers 



As = -^u' 



^--4(,;.„>P.«,).^- 



2A4 



■N' 



6N' 



5N' 



(6.34) 



dk{Juu'P^,N,]+^^ 



3.35) 



Substituting these back into (6.33), we end with the field 
equation for the director field: 






dkiJuu'^P'.Ni 



5V 



3.36) 



Here iTi^ is the projector defined in (5.20). Note that we 
cannot yet put equal zero the left-hand side of (6.36); it 
will be evaluated later from the variation of the electro- 
magnetic part (5.30) of the total Lagrangian. 

Contracting (6.32) with u', we find another Lagrange 
multiplier: 



2An 




1 



tu'N, 




3.37) 



Here again we cannot put equal zero the first two terms 
on the right-hand side of (6.37); they also should be in- 
serted from the variation of (5.30). 

Finally, let us notice that the material Lagrangian "on- 
shell" (i.e., after making use of the field equations) reads: 



C 



p + v- 



SC' 



5v 



-V. 



(6.38) 



VII. CANONICAL NOETHER CURRENT 
TENSOR FOR MATTER 

The general form of the canonical energy-momentum 
tensor for the matter part is 



dC' 



ddj^"^ 



a-*^ - 5,:> r 



(7.1) 



Here "^^ are all dynamical variables. The matter La- 
grangian (6.1) depends on the derivatives of A^' as well 
as on Ai, s, and X. The derivatives with respect to the 
velocity of the director field are computed from (6.8) and 
(6.12): 



dC 



dC^ 



^ dCP ^ _ k pj jy _ dV 

ddkN' ^ ddkN' ^ ddkN^ ~ """ ' ' ddkN^' 

(7.2) 
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On the other hand, from the constraint Lagrangian we 
find 






ddiS 



A2U\ 



dd,X 



A3U^ 



(7.3) 



Consequently, we have 



dV 



a■iv^ 



(7.4) 



and 



dd,A, 



a, Ai + TT7^9,s + ■7^TrT7djX = u' {-vdjKi 



dd,s ^^ 



dd,X 



'A2djS + A3djX). 

(7.5) 



In order to compute the right-hand side of (7.5), we can 
use (6.27). The latter equation yields, with the help of 
(6.37) and (6.34): 



77 ■ 



SC 



J 



-P. 



+ P + P+-UUJ 

ov 2 

1 . \ . A/"'" BV 

Jv \ PrdkN' - - u^'Nk ]Nn + ^^ + jrT 



-N. 



^C2 



6C^ 



dn (jvu^' 



PkNi + 



5V 

(7.6) 



Substituting (6.38), (7.4), (7.5), and (7.6) into (7.1), we 
find explicitly the canonical energy-momentum tensor of 
the liquid crystal: 



Here the Frank deformation stress tensor is 



F . 



dV 



diN^ + (Sf V, 



the effective pressure p°^ reads 



„cff 



p + V- 



SC 



5v 



(7.7) 



(7.^ 



(7.9) 



and the relativistic 4-momentum density of the fluid is 



1 / 1 T ' 



~P' 






-Nk u^Ni + — ■ 






dw' 



dn Jvu'^Pk^Nl + 



5V 



5N'' 



(7.10) 



As usual, we have to compute the variational derivatives 
of the matter Lagrangian with respect to v^u^^N^ from 
the electromagnetic part of the total Lagrangian (5.30). 



A. Canonical spin of liquid crystal 



Denoting all the fields in the system by the symbol 
^'*, that carries a "multi-index" A, with {pij)'^B as the 
generators of the Lorentz algebra for these fields, the spin 
of the system is defined by 



S,'-J^ip^,)i^''. 



d{dk^^) 



(7.11) 



Since the matter Lagrangian does not depend on the 
derivatives of u', the spin of the system is straightfor- 
wardly computed with the help of (7.2): 



S^■i' = N,, 



dC 



ddkNi\ 



= -Jvu^N^iP'j^Ni -N\ 



dV 



ddkNi] ■ 
(7.12) 



Using the results of the Appendix A, we can write down 
the variational derivatives for the Frank potential (6.4). 
They read; 



dV 
dV 



Ko 



{^^N, - djN,){N^ - u'^d'Nk) 



- -{K2 - K3)iNPdpN,){N''dgN'')uk, (7.13) 



= KiidkN'')di + K2 PlPlid'^Ni - diN'') 



+ (K2-K3)N^P^^{NPdpNk) 
SV 



(7.14) 



SN' " ^' 



Ki{dkN'')6l + K2 P^Plid'^Ni - diN'') 

+{K2-K3)N^Pt{NPdpNk)] 
+ {K2 - K3mNP){N''dkNg)P^. (7.15) 

As a consequence, the Frank stress tensor (7.8) is 

- K2{d,Nk)p^p^{dpm - dmp) 

- (K2 - K3){d,Nk)N^Pl'{NPdpN') 






(7.16) 



and the spin density tensor of matter (7.12) reads 

^./ - - Jvu'^N^.P',^ Ni - KiN[, 6^^ diN' 

- K2P''''%P'j]{dnNi - diN^) 

- (A'2 - if3)^^A^[,^',]A^"a„7Vz. (7.17) 



B. Balance equation of the angular momentum of 
matter 

Let us check that the angular momentum balance 
equation is satisfied for the open material system under 
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consideration. This balance equation follows from the 
Noether theorem and it reads [2Gl: 



S[y] + dk Sij = - 






A,T,B 



S-qjA 



(7.18) 



In our system, we have two vector fields, w* and iV*, and 
the corresponding Lorentz generators are (p'j)S = SqS^^ 
As a result, the right-hand side of (7.18) explicitly reads 



SC 



iPv)>' 



SC 



6uP ^"^'"^ 



SNP ^^'^'"^ 



ip^,r,N'' 



SC 



6C' 



A^,i 



(7.19) 

From (7.7) we derive the first term on the left-hand 
side of the balance equation (7.18): 



Denote 






(7.20) 



(7.21) 



then the field equation for the director (6.36) is recast 
into 

^.'■*fe-0. (7.22) 

Let us consider the last term in (7.20). By making use 
of (7.10), we have 

SC' . dV 



r[, uj] 



r- Uil H r- Mil 



^£™ dV 

-iV[,$,]+7V[,7r,]'=$,. (7.23) 

We used the definition of the projector (5.20) that yields 



tt/' -S'' - NjN'' 



(7.24) 



to transform the two last terms on the first line in (7.23). 
In view of the field equation (7.22) the last term in 
(7.23) vanishes, and we find 

dV 
dul- 



6C" 



dV 



Nr 



5ui^ ^' 6Ni^ '^ ' du^^ ' 



SV 



A^oi 



JzyA^[,P%]7Vfc 






(7.25) 



Recalling that Ni = u^dnNi, we see that the last line 
reduces to the total divergence 



^n 1 t- 



7V[,a„ ( J:/M"Pj]'=iVfc = dn {jiyu''N^,Pj]''Nk 



. k aV. 



(7.26) 



After these preparations, we can find the left-hand side 
of the balance equation (7.18): 



S[ij] + dk Si 



SC" 



5C" 



■N, 



5ui^ "■" ' 5Ni^ ''^' ' ^'^^'l 
SV ( dV 



dV 

. (7.27) 



The last four terms in (7.27) disappear in view of the 
Lorentz invariance of the Frank potential V that results, 
via the Noether theorem, in the following identity (see 
appendix A): 



F dV 



SV f dV 



= 0. 

(7.28) 



Finally, using the identity (7.28) in (7.27), the angular 
momentum balance equation for matter can be recast 
into 



5C" 



Sfiji + dk Si] ' — -T- 

^ •'^ Su 



(5£'° 

sm ^'1 



(7.29) 



Accordingly, by comparing (7.29) with (7.19) we indeed 
verify the identity (7.18) of the angular momentum for 
the open material system. 



VIII. CANONICAL ENERGY-MOMENTUM 

TENSOR FOR THE ELECTROMAGNETIC 

FIELD INSIDE THE MEDIUM 



The canonical energy-momentum tensor for the elec- 
tromagnetic field is obtained with the usual definition 
[26] applied to the electromagnetic Lagrangian C^™ in 
(5.1). If we consider the 4-potential covcctor Ai as the 

fundamental electrodynamical variable, then the electro- 
cm . 
magnetic canonical tensor S i^ reads 



em . dC^™ 



S,^ := 



d]Ak -S]C 



dd.Ak 

E^+i7^^(9feA0 



(8.1) 
(8.2) 



M 



where Ei"* is the usual Minkowski tensor for the electro- 
magnetic field in matter, given by 



M . ., 

S.^ := -F.kB^'^ 



-A'FkiH^'. 



(8.3) 
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Now, using the expression (5.29) in (8.2) and (8.3), we 
get explicitly 



M ■ 1.1 IX 



4 



1 



Moc^ 



(e± - A^l' - A/i-i) [-F^Vi^,;w' 



Mo 






^iY^F^feF'-^'ATz] 



Moc^ 



(Ae + Afi-^) {FpgNPt 



-5l{FuN^u') 



-u^F.^N-^ + iV-'F„M"] , 



(8.4) 



which is the explicit expression for the Minkowski tensor 
of the field inside the liquid crystal. 



electromagnetic Lagrangian inside matter £*"" describes 
an open system. 

On the other hand, computing the variations of £'^'" 
in (5.30) with respect to the material variables yields. 



6L° 



1 



1 



Hoc 
(5£™ 1 



(Ae + A^-') (FpgNPui) F,fciV^ (8., 



SN' Mo 



A^i-^F,kF'''Nl 



+ J (Ae + A/i-i) (FpgNPu'^) F.^u'^. (8.9) 

Comparing (8.8) and (8.9) with (8.7), we immediately 
verify the correct balance equation for the electromag- 
netic angular momentum part of the system. 



cm j^ J/:™> SC 

S [ij] + Ok S ij — „ r- Ui 



<5u[» "^' ' <57V[' ^^- 



(8.10) 



A. Balance equation for the angular momentum of 
the electromagnetic field 



This is in perfect agreement with the general Noether 
identity (7.18). 



From the general definition [2C], the spin density S ij'' 
of the electromagnetic part of the system is given by, 



dC'' 



didkAr^ 



-ip,j)jAi^H\,Aj]^0. (8.5) 



IX. TOTAL CANONICAL 
ENERGY-MOMENTUM TENSOR 

The complete system of material medium plus electro- 
magnetic field is described by the total Lagrangian 



Now we are in position to evaluate the angular momen- 
tum balance equation for the electromagnetic part of the 
system, which has the same form as the one for the mat- 
ter part (7.18). Taking the antisymmetric part of (8.2) 
and using the expression (8.5), together with the Maxwell 
equations without sources, djH"^^ = 0, we see that the 
left-hand side of the identity for electromagnetic angular 
momentum is simply given by the antisymmetric part of 
the Minkowski tensor: 



cm , M 



S [ij\ + dk S ij — S 



(8.6) 



where 

M 

Sri,-] : 



U[jF,]k 



Moc 



^ {e± - Ml' - A/i-') 



F^'ui 



MoC 



{Ae + Afi-^) (FpgNPu'i) N' 



N[jF,]k 



1 

Mo 



-Ah^'F^'Ni 



+ 2 {Ae + A/i-i) (FpgNPu'^) u^ 



Moc 



(8.7) 



C := £'" + C 



in I ncn\ 



(9.1) 



As a result, the total canonical energy-momentum tensor 
of the closed system is given by 



in . cm . 



(9.2) 



with the electromagnetic part Si-' given in (8.2) and 

111 
(8.4) and the material part Yn^ given in (7.7). 

In order to find an explicit expression for the total 

canonical energy- momentum tensor Ei^ , we first need to 

evaluate the variations SC'' /5v,5C^ /5u' and 6C^/5N' 

and then insert them in (7.9) and (7.10). For this aim, we 

take into account the equations of motion of the material 

variables 



5C^ 

5cr_ 



SC' 



5v 



0, 


(9.3) 


0, 


(9.4) 


0, 


(9.5) 



We see that the electromagnetic canonical energy- 
momentum tensor as well as the Minkowski tensor are 
not symmetric. However, it is not surprising that the 
right-hand side of (8.7) is not equal to zero since the 



from where we clearly see that the variations of the mat- 
ter Lagrangian are exactly the negative of the variations 
of the electromagnetic Lagrangian, which we have al- 
ready explicitly computed in (8.8) and (8.9). In addition. 
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the variation of the electroniagnetic Lagrangian with re- 
spect to the particle number density i/, exphcitly yields 



SL'' 



5v 



1 ^M± g2 



2 \ dv MomI dv 

2 \ ou fj,o ov 

(9.6) 

where we defined the 4- vectors electric field £i and mag- 
netic field B^ in (5.31) and (5.32), respectively. 

The variation (9.6) enters in the expression of the 
"effective" pressure (7.9), which include the terms de- 
scribing the electrostriction and magnetostriction effects. 
Then, replacing the negative of the three variations (8.8), 
(8.9), and (9.6) into (7.7), (7.9), and (7.10), we explicitly 
obtain 

+ -^{Ae + A^i-^)u^P.i^FMN\Fj,,N^u'i), 

fXQC 



where we denoted 



— u,{p^JvLo''/2)+P,'' 



0, := dk (^Ji^u'^P^N,^ 





\ Jv <r 






' 1 


+ F,'^ 




du^ 


' 








(9.8) 


SV 








SN'' 






(9. 


9) 



Since we already checked that the angular momentum 
balance equations are fulfilled both for the matter and 
electromagnetic parts of the total closed system in (7.29) 
and (8.10), respectively, it is obvious that if we add both 
equations, then the angular balance equation for the total 
system will be also valid. 



<9fc5„'^ = 0. 



where 



^ii •'~' Ji 



O i 



(9.11) 



(9.12) 



It is worthwhile to notice that the right-hand side of 
(9.11) vanishes since the total system is closed; however, 
the total energy-momentum tensor (9.10) is not symmet- 
ric, since the spin density of the system (9.12) is nontriv- 
ial. 



X. FULLY EXPLICIT ENERGY-MOMENTUM 
CONSERVATION LAW 



The total system under consideration, composed of a 
rclativistic liquid crystal plus electromagnetic field, is a 
closed system. There are no external fields present like 
J*xt and therefore the total canonical energy-momentum 
tensor (9.10) of the system is conserved: 



djll^ = 0. 



(10.1) 



The definitions oiTi-' andp'' are given in (7.8) and (7.9), 
respectively. Finally, if we consider (8.2), (8.4), and (9.7), 
we can write the total canonical tensor (9.2) explicitly, 
which reads 



f (mI^ + Am-^) 



1 

4^ 



F'^'F^k + ^SiF'^'FM 



1 



MoC 



(e± - /il' - A/i-i) [- F^'^ukFuu' 



/"O 



-Afi-^ 



+ ( ^Si - ^u,u^ ] {Fuu'f 



F^'^NkFuN' + ^SfiFkiN'f 
+ NJF,kF'''Ni + ^N.u^u'^FkiF'^'N^ 



+ 



Hoc' 



{Ae + A/i-i) {FpgNPu'^) [iV^F^u" 



+ ( - \si + ^u,v? ) {FuN^u') 



H'^KdkA,). 



(9.10) 



Notice that in the expressions (7.12), (8.5), and (9.10)- 
(9.12), we have obtained the total energy- momentum ten- 
sor of the system and the identity of total angular mo- 
mentum of the system, but without using a specific ex- 
pression for the potential function V = V(A^*,9j7V%m*). 
As a consequence, all the latter expressions are valid 
for an anisotropic uniaxial diamagnetic and dielectric 
medium, with any internal dynamics for the 4-director 
field N"^ . By considering the expression (6.4) for the 
Frank potential V, one can derive explicit expressions of 
the derivatives 5V/5u\ dV/d{djN^), 5V/SN' and of the 

F . m 

tensors Tj-', Sij ■ This is done in detail in Appendix A. 
In order to obtain an explicit expression of this con- 
servation law, we can insert (9.8), (9.9), (7.13), (7.15), 
and (7.16) into (9.10). Due to the macroscopic Maxwell 
equations (5.2), the gauge non-invariant term in (9.10) 
vanishes and we finally obtain 



d,j:^ + 9jS,^ = 0. 



(10.2) 



Here the non-symmetric energy-momentum tensor Y,i-' 
only depends on the material variables, as if the rela- 
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tivistic liquid crystal were in isolation, 



where we defined 



1 



■ UiU- 



p-^Jiycj'^ojk] -P^'p''^ 



h, := - , . 



5V dL" 



dNi 



(11.2) 






K2 {dkNi - diNk){N' - vP&Np) 



-JvNk u'Ni + {K2 - K3)iNPdpNk)iNWgN')ui 



as the total 4~i^olecular field, since its spatial compo- 
nents reduce, in the non-relativistic limit, to the standard 
"molecular field" [1, 16]. In Sec. XII we study this limit 
-n more detail. In order to better interpret the dynamics 



Ki{dkN''){d^N^) - K2{d^Nk)P^Pq{dPN'^ - d'^N^^i N\ we can contract (11.1) with eP«*iV„ make use of 



(i^2 - K3MNk)N^Pt{NPdpN') 



Ki5^ {dkN'^y 



-20i- 



^kln 



'NkdiN„ 



KsS^ (ep„,7V™e'='"9,7V„ 



Jiyuku'di P/W^ +XiM'^afe(a™iV") 



:N,ul 



+ {K3- K2)uHdkNP){N'diNg)Pp^] , (10.3) 

except for some electromagnetic terms inside the effective 
pressure p*^*^, which describe possible clcctrostriction and 
magnetostriction effects: 



„cff 



^1 f de ^, 



1 di-L± 



B' 



i\2 



^ dv 
1 5A/1-1 



A^o 



di 



(10.4) 



The second term in (10.2) corresponds to the 4- 

divergence of a non-symmetric "field" energy-momentum 

f ■ . 

tensor Sj-', which contains electromagnetic field terms 

coupled to the material variables: 
1 , 



f . 



Mo 



1 



A/i" 



iei_ - At_ 






s^^f'^^f, 



kl 



'-F^^UkFuu' 



S' 



1 



l\2 



{Fkiu') 



Mo 



-A/i" 



- F'^^NkFaN' + -Sl{FkiN'f 



+N^F,kF^'Ni + -N.u^u'^FkiF'' 



'AT 



+ 



Moc^ 



(Ae + A^i-') {FpqNPu'i) [WF.nu'^ 



- T^H + -u,u> {FkiN'^u') 



(10.5) 



XI 



RELATIVISTIC DIRECTOR DYNAMICS 

Let us analyze the equations of motion for the director 



N^. Rewriting equations (6.36) or (7.22), we have 



d. 



Jv u^ P'j Ni 



= 0, 



(11.1) 



(A6) and the continuity equation (6.19), to obtain 



Jv-k"- 



iW^ =-€'''' Njhk. 



(11.3) 



From (11.3) we see that the 4-molecular field hi is re- 
sponsible for the "torques" and changes in the 4-director 
A^* of the liquid crystal. We can identify two contribu- 

F cm 

tions to the 4-molecular field hi :^ hi + hi- One is the 
Frank deformation 4-molecular field 



F 

hi '■= — 



SV 

In* 



(11.4) 



= d, [Ki{dkN'')S^ + K2 Pk'PHd'^N, - diN'') 
^{K^-K2)N^P^\N'P^pNk)\ 
+ (^3 - K2MNP){N^dkNg)Pp\ (11.5) 

which describes the changes in N"^ caused by the defor- 
mations of the liquid crystal itself, and the other contri- 
bution is the electromagnetic 4-molecular field 



h^ := 



dC" 



— A/i-i 
Mo 
1 



(11.6) 



FrkP'^'m 



+ 2 (Ae + A^-i) {Fp.NPu'^) F,ku\ (11.7) 

Moc^ 

which describes the infiuence of the dynamical electro- 
magnetic field on the orientation of the 4-director. 

With the knowledge of (10.2), together with Eqs. 
(11.3) for TV', the continuity equation di{yu'^) ~ in 
(6.19) for u% V and the Maxwell equations (5.2) for Fij, 
we can completely determine the dynamics and evolution 
of this system, composed of the relativistic liquid crystal 
with anisotropic optical properties interacting with the 
electromagnetic field. 



XII. NON-RELATIVISTIC LIMIT 

Let us study the dynamics of the liquid crystal in the 
non-relativistic limit, when the motion of the fiuid is such 
that |v| ^ c. In particular, this approximation can be 
applied when the liquid crystal is at rest in the laboratory 
frame. In the non-relativistic limit, we expect consistency 
with the earlier results [1, 8, 16, 17], but first we need 
some technical preparations. The 4-velocity reads 



u' ^l{l,v), 



(12.1) 
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and therefore we have for the components of the projec- 
tor: 



p b _ c b 



Po" 



l'v\ 



pQ _ 2^ 



Po°=7' 



(12.2) 
(12.3) 



Hereafter, the 3-dimensional indices are raised and fow- 
ered by the EucHdean metric; in particular, v = v"', 
Va = Sabv'', w^ = v°'Va = Sabv'^v^, ctc. For the skew- 
symmetric tensor (4.5) we find expUcitly 



^abc ^ ^^abc^ 



^abc 'J^abc-! 



^Qab 






eOah = l<^abcV 



(12.4) 
(12.5) 



Taking into account the orthogonality conditions (4.2)- 
(4.4), the 4-director TV' reads in components 



iv^=('(!i:;i),„ 



(12.6) 



Notice that the 3-vector n, with Cartesian components 
n", recovers its normalization n^ = dabn""^'' = 1 only 
in the non-relativistic limit. In general, it satisfies n^ = 

1 + {v -nY l(?. 



A. Non-relativistic director dynamics to zeroth 
order in vjc 

First we note that iVjCJ-' = in (11.3) and therefore we 
obtain 



taking the non-relativistic limit to zeroth order in vjc of 
(11.7), which explicitly reads 

cm , ^ 

ho :=^ eoAeiE ■ n)E - fi^ Af^ {B ■ n)B 

+ fi^^Afi-'^B^n. (12.11) 

The result (12.11) looks slightly different from the usual 
electromagnetic molecular field of the non-relativistic 
models [1, 8, 17], since we use the fields E and B as the 
independent fields and not E and H . One can, however, 
easily recover the same formulas in a different disguise 
by making use of the inverse of the constitutive relations 
(2.2). Notice that the last term in (12.11) is proportional 
to the director n and therefore it does not contribute to 
the torque when replaced in the cross product of (12.8). 
Therefore, we can ignore this last term and redefine the 
electromagnetic molecular field to zeroth order in v/c, as 

c 

a sum of two terms, the electric molecular field ho and 

m 

the magnetic molecular field ho , given by 



ho = eoAeiE ■ n)E, 
ho= Ati-\B-n)B. 



(12.12) 
(12.13) 



Collecting together all the results of this section, the 
dynamics of the director n in the nonrelativistic limit, to 
zeroth order in v/c, is described by 



Jvny. 






n X h„ +n X h': +n X h" 



(12.14) 



where the molecular field has independent contributions 
from the Frank deformations and the interactions with 
electric and the magnetic fields. 



JiyPi'cu^ = -e'^'^Nih 



jik- 



(12.7) 



Then, inserting (12.1)-(12.6) into (12.7), we find the non- 
relativistic equation for n" to zeroth order in v/c 



Jv6jo = n X ho- 



Here, 



ujo ^ nx 






(12.. 



(12.9) 



is the angular acceleration of a liquid crystal fluid element 
to zeroth order in v/c. The right-hand side of (12.8) is 

F cm 

determined by the molecular field, ho := ho + ho- The 
fiuid part 

L° = {Ki - K2){dadbn'') + A'a^^abV^n" 
+ (Ks - K2)Sad dbin'n-dcn'') 
- {K3 - K2)SM{dan''){n'-'dcn''), (12.10) 

is the zeroth order Frank deformation molecular field (cf. 
the non-relativistic equation (2.147) of [8]). The electro- 
magnetic part of the molecular field can be computed by 



B. Non-relativistic solutions 

Let us assume that the electromagnetic field vanishes, 
so that hl = h'^ = 0. 

Stewart in Sec. 2.5 of [s] describes some exact solutions 
of the equations of motion. One is an obvious constant 
director solution; the other is the static spherical solution 
and the twist solution [with the local coordinates x = 



(cos 6*, sin 0,0) . 



9 = cix^ 



C2- 



(12.15) 
(12.16) 



Here we notice that the static twist solution (12.16) can 
be generalized to a dynamical "plane wave" solution: 

n'' = (cose, sin 8,0), O ^ cot + cix^ + c-i, (12.17) 

where co,ci,C2 are arbitrary constants. In this case, we 
explicitly have. 






h' 



K2{cifn, 



(12.18) 
(12.19) 
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and thus (12.14) is satisfied. In the Uterature, the wave 
solutions of the full nonlinear equations of motion has 
attracted some attention (see [16, 17, 56]). 



XIII. SUMMARY AND DISCUSSION: 
ABRAHAM-MINKOWSKI CONTROVERSY 

In this paper, we have constructed a complete relativis- 
tic Lagrangian theory of a ncmatic liquid crystal. Our re- 
sults provide a consistent relativistic model for a medium 
with anisotropic optical properties in interaction with the 
electromagnetic field. In particular, in such a framework 
one can study the problem of the proper description of 
the energy and momentum of light in anisotropic media. 
This should shed light on the long standing Abraham- 
Minkowski controversy, traditionally discussed only for 
isotropic media. 

We have generalized the earlier non-relativistic model 
[16, 17] and the variational model of an ideal relativistic 
fluid [2()] and derived a complete theory of the nematic 
liquid crystal medium and its interaction with the elec- 
tromagnetic field. We have derived the nonlinear equa- 
tions of motion for the liquid crystal fiuid, and explicitly 
verified the total energy, momentum, and angular mo- 
mentum balance laws, which arise as consequences of the 
Noether theorem from the invariance of the (field plus 
matter) system under spacctime translations and Lorcntz 
transformations, respectively. In this work the general 
formalism is presented in full detail. We will analyze the 
solutions and applications separately. The analysis of the 
properties of electromagnetic waves in the moving liquid 
crystal will be also considered elsewhere. 

As we have seen, liquid crystals are an interesting ex- 
ample of continuous media with microstructure. The 
"internal" degrees of freedom of such a medium is repre- 
sented by the director vector field iV* assigned to every 
material point of the fluid. This field gives rise to a non- 
trivial spin of the medium (7.17). As a consequence, the 
total energy-momentum tensor (9.2) of the closed system 
of the medium plus the electromagnetic field is not sym- 
metric. This asymmetry is crucial for the validity of the 
Noether identities. 

In this paper, we have neglected dissipation. The gen- 
eral formalism developed here is applicable to any mov- 
ing medium with uniaxial anisotropic properties. In par- 
ticular, one should have in mind possible astrophysical 
applications [-54], where our model provides an explicit 
dynamical mechanism for the description of a physical 
medium with uniaxial anisotropy. As concerns the liq- 
uid crystals which belong to the class of moving uniax- 
ial media, they are characterized by a nontrivial dissi- 
pation. In this sense, our theory has limited applica- 
bility and it should be considered only as a first step 
in constructing a full realistic physical model. The rel- 
ativistic mechanics of dissipative fluids has a long and 
controversial history, with the flrst attempts going back 
to Eckart ["i7] and Landau-Lifshitz [iS]. Later an es- 



sential improvement was achieved in the works of Israel 
and Stewart [59, 60]. However, these models used a phe- 
nomenological approach with numerous ad hoc assump- 
tions, and it seems more appropriate to use the approach 
of Carter [()1-G3], which generalizes the variational prin- 
ciple by replacing the Lagrangian with the so-called mas- 
ter function that systematically takes into account the 
irreversible viscosity and thermal effects. The develop- 
ment of the Carter type approach for dissipative media 
with microstructure will constitute the next step for the 
construction of the variational theory of liquid crystals. 
An important feature is that the Noether type identities 
play a central role in Carter's approach. 

Our current results contribute to the discussion of the 
energy and momentum problem of the electromagnetic 
field in a medium. In particular, our analysis clearly 
demonstrates that for the case of an uniaxial anisotropic 
medium it is the total energy-momentum tensor of the 
coupled system (matter plus field) that is important for 
the understanding of the balance of the momentum and 
angular momentum. For the case of an isotropic medium, 
this was pointed out more that 40 years ago by Penfield 
and Haus [28, 29] and more recently in [26, 27, 30]. 

Notwithstanding this general fact of the importance of 
the total energy-momentum, one can split it into a "mat- 
ter" and a "field" part in many different ways. Specifi- 
cally, by expressing (9.10) in terms oi N^, u^ , Fij, permit- 
tivities and impermeabilities, we have shown that there 
exists a split of the form (10.3) plus (10.5). The purely 
matter part (10.3) does not depend explicitly on the field 
strength Fij, except for the terms present in effective 
pressure. Therefore, it could be identified with what is 
sometimes called a "kinetic" energy-momentum [31, 32]. 
The corresponding field part (10.5) has a form which 
is more complicated than the usual Abraham tensor for 
isotropic media [2(i], since the former involves not only 
the field Fij, the 4- velocity m*, the isotropic permittivity 
e and the permeability ^, but also the 4-director iV', its 
derivatives, and the anisotropics Ae, A^~^. Despite the 
fact that its structure is different from that of Abraham, 
this tensor plays a role analogous to that of the tradi- 
tional Abraham tensor and moreover reduces to it in the 
case Ae = A/i~^ = 0. 

However, the tensor (10.5) is not symmetric. It was 
pointed out in [2()] that, for simple media where the 4- 
velocity u' is the only non-scalar field contained in the 
constitutive relation, the total energy-momentum tensor 
turns out to be the sum of a kinetic term plus the Abra- 
ham tensor. It was not clear whether the same is the 
case for anisotropic media. Our results now clarify this 
point. Another interesting problem would be to find a 
consistent definition of some kind of "generalized Abra- 
ham tensor" , possessing the same properties of the Abra- 
ham tensor in the isotropic case, but valid in any type 
of medium. This would certainly improve our under- 
standing of the Abraham-Minkowski controversy, which 
has been restricted only to simple media. Results in this 
direction will be analyzed in a forthcoming publication. 
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Appendix A: Explicit calculation of the Frank 
potential terms 

In this appendix the exphcit derivatives of the Frank 
potential are computed and an exphcit verification of the 
Noether identity (7.28) is given. For the computation of 
the derivatives of the Frank potential (6.4), the following 
identities will be useful: 



p. ., fPi^ — pp pi. ^ pp.pg . 



(Al) 
(A2) 

where P^ j and e^-fc are defined in (6.7) and (4.5), respec- 
tively. Then, using (Al)-(A3), we can prove these other 
useful relations: 

e^jkN^e'^'^diN,, = P^N^diN^, (A4) 

pi^.^VJ ^e*J"a;jiV„, (A5) 

2iV[fce,],,^'iV^' = - e,,„w" = ^N^^P^^ N,. (A6) 

We now are in condition to compute the contribution 
of the Frank potential to the dynamics of the fluid. For 
this, it is convenient to separate the potential (6.4) into 
three pieces: 



V:=Vi+V2+V3. 



The first elastic constant 



(A7) 



For the splay deformation elastic potential. 



Vi = \K,{d.,N^)\ 



we find 

dVi dV 



= 0, 



ddjN 



-=K,{duN^)5l 



9Vi .. ,. .,, 



(A8) 

(A9) 
(AlO) 



Then, using (A9)-(A11) in (7.8), we obtain 



Fl 



T^' = - K^idkN^mN^) + SlVi 



(A12) 



Fl 



Tfo-] = - A-i(afciV*^)(%iV,]). (A13) 

Finally, let us calculate explicitly the terms: 

^^N,^^-K,N^d,^{duN\ (A14) 
dVi 



'dk A^[,; 



ddkN^i 



-^~ KiidkN'')idyN,^) 

- A'iA^[,a,](5fe7V^). (A15) 

Substituting all this into (7.28), we verify the Noether 
identity for Vi. 



2. The second elastic constant 

The second term (twist deformation) in the Frank po- 
tential is given by 



V2 = lK2{e''''md,Nk)\ 



(A16) 



Using the identities (Al) and (A3), we can rewrite this 
as 



1 



V2 = -K2 Kirf {d,N, - d,N{){d''N' - d'N^) (A17) 



= -A'2 
2 



PlPl{d,N,~d,N,)d''N'' 
+P]{NPd^N,){md,N^)\ . 



(A18) 



The derivatives of the potential V2 with respect to the 
director and its derivatives are straightforwardly calcu- 
lated: 

ay; 

^ = K2{d.Nn{N''dkN,)P^^ 

= K2{d,N'^){NPdpNk) 



^{diNP)up{N''dqN'')uk, 



(A19) 



dV2 
ddjN' 



^K, 



PlPlid^Ni - diN") + N^PtiNPdpNk) 



= K2 {d^N, - d^N^ + N^NPdpN,) 

_ :^ U(N, - u''d,Nk) - u,{N^ - u'^d'Nk) 

+N^Ui{NPdpN'')uu\ . (A20) 

In addition, the derivative with respect to the velocity 
reads 



dV2 



K2 

r2 



+{NPdpNi){N''dqN'')uk] . (A21) 



Substituting (A20) into the definition (7.8), we obtain 
the contribution of K2 to the energy-momentum tensor 



F2 



T^ = 51 V2 - K2{d,Nk) [p^p^id^m - d^P) 

+N^Pl'iNPdpN^)] (A22) 

= Sf V2 - K2{diNk) [d^N'' - d^N^ + N^NPdpN'')] 



& L 

+ N'\NPdpN^)u^uA . 



(A23) 



In order to check the Noether identity (7.28) for the 
V2 term, we first notice that it can be identically recast 
into 



F2 



5V2 



^M + a^^j] 



5V2 9V2 

^^iV,.]-(a,iV[,)^^^^ 



= 0. 

(A24) 
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Using (A20), we derive an intermediate result: 
dV2 



{dkN,) 



ddkN' 



= K2[{d''Nj){dkN,^d,Nk) 



n2 



-{NPdpNj){NWqN,)] 



~u^{d^Nj){Nk-u^dkNi) 
+u,u\N'PdpNj){mdqNk)\ . (A25) 

It is straightforward to find the antisymmetric objects 
using (A19), (A21), (A23), and (A25): 

F2 , K-) ■ , 

+ K2N[,{d,]N''){NPdpNk) 



K 



- -^Ui,idj]Nk){N'-u'd''Ni) 



K- 



-^N[,{d^]N'')uk{mdgNi)u\ (A26) 



. K2 



■ U[,NP{d\p\N^^)iNW,Nk)u\ (A27) 






-TVj] = - K2N[,id,^N^){NPdpNk) 



K- 



+ -^ N[,{d^]N'')uk{mdgNi)u\ (A28) 



and 



^^'^^'^ W^\ ^ i^2(a'=%-)9,]iVfe - ^Ny^{d,^N,)u'' 



(A29) 

Substituting all this into (A24), we verify the Noether 
identity (7.28) for the second term V2. 



3. The third elastic constant 

Analogously, we consider the third term (bend defor- 
mation) in the Frank potential, 



V3 = - -K3 (e,,fc7V^e'='"a,iV„) 



(A30) 



K3P.^{NPdpNj){N''dgN'), (A31) 



where we have used (A4) to simplify the potential. It 
is worthwhile to notice that this quadratic invariant has 



the same form as the last term in the K2 potential (A18). 
As a result, the corresponding derivatives of the potential 
(A30) with respect to its arguments can be conveniently 
extracted from the formulas (A19), (A20), and (A21). 
These derivatives are given by 



dV 



- K3 {d,N''){NPdpNk) 



(A32) 



^{^^NP)up{NWgN'')uk, (A33) 



ddjNi 



- KsN^PtiNPdpNk) 
= - K3 N^NPdpN, 



K. 



N^u,{NPdpN'')uk, 



dVs K3 



iNPdpN,)imdgN'')uk. 



The stress tensor (7.8) for V3 reads 



(A34) 

(A35) 
(A36) 



F3 



T^' = <5^' V3 + K3[d,N,,)WPf{NPdpN') (A37) 



^, V3 + K3{d,Nk)W{NPdpN'') 
- ^{d^N,)N^NPdpN')u''ui. 



(A38) 



In order to check the Noether identity for this last part of 
the potential V, we observe that Eqs. (A33)-(A38) yield 



F3 



T[.,] = - K3N^(d,^N''){NPdpN^) 



K: 



+ ^ N[,{dj]N^)uk{N'^dgNiy 



|§«.] = - ^U[,NPid^p^N^^){N''d,N,)u'^ 



^^N,^=K,Ni,{d,^N''){NPdpN,) 



(A39) 

)u^ 

(A40) 



idkNu) 



dVs 



K, 



K. 



N[,{dj]N'')uk{mdgNi)u' 



(A41) 



ddkN^] 



u^,NP{dip\N^^)iN'^dgNk)u'' 



(A42) 



Substituting all this into (A24), we verify the Noether 
identity (7.28) for the third term V3. 
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